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Abstract 

We consider a quantum system of non-interacting fermions at temperature 
T, in the framework of linear response theory. We show that semiclassical 
theory is an appropriate framework to describe some of their thermodynamic 
properties, in particular through asymptotic expansions in h (Planck con¬ 
stant) of the dynamical susceptibilities. We show how the closed orbits of the 
classical motion in phase space manifest themselves in these expansions, in 
the regime where T is of the order of Ti. 
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1 Introduction 


Consider a system of non-interacting fermions confined by an external poten¬ 
tial and in contact with an exterior reservoir at temperature T. Assume that a 
time-varying external perturbation drives the system out to, but near of, its equi¬ 
librium state. The response of this quantum system to an external time-dependent 
perturbation is a subject of high physical interest, which can be investigated experi¬ 
mentally, in particular the so-called “dynamical susceptibility”. A complete rigorous 
analysis of this problem is still lacking, although recent progress is being made in 
the understanding of non-equilibrium statistical mechanism, and its link with the 
underlying chaotic dynamics 00 0. 


A semi-empirical route which has been proposed (see classical textbooks |lj, [15 


consists, for small perturbation, of investigating the response function “to first order 
of the perturbation”, i.e. the so-called “linear response theory”. This semi-empirical 
route has been given a firmer foundation (see the book by Bratclli and Robinson 
0) where a link with the KMS condition is established. (See also recent progress 
in ®). 

In this paper we rederive the first order response function for the quantum 
fermionic system under study, i.e. the so-called “generalized Kubo formula” (see 
also 0) and investigate semiclassical expansions of it, assuming suitable “chaoticity 
assumptions” on the one-body underlying classical dynamics. These semiclassical 
expansions are developed in a similar spirit as previous studies on the “semiclassi¬ 
cal magnetic response for non-interacting electrons” ]!]. [|, [10], [L3], [L6], [17|] i.e. we 

exhibit a low temperature regime where the closed classical orbits of one-particle 
motion manifest themselves as oscillating corrections to the response function. 

Section 2 contains , for pedagogical purposes the basic framework of so-called 
“second quantization” in which the physical system under consideration can be 
studied and its thermodynamical properties mathematically investigated. Section 
3 presents the so-called “linear-response theory”, and the dynamical susceptibility 
that will be studied in the semiclassical framework. Section 4 presents and derives 
the main results of this paper: a rigorous semiclassical expansion of the dynamical 
susceptibility under suitable assumptions on the physical system. 
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2 The physical model 

Consider a system of non-interacting fermions, living in IR n , subject to a one- 
body Hamiltonian H which is the Weyl quantization of a classical Hamiltonian 
H(q,p) of the form 

2 

H(q,p) = ^- + V(q) (2.1) 

2m 

with V G C°°(R n ) such that the following confining assumption holds: 

Assumption 1: 


V(q) > Cq(1 + q 2 ) s/2 s,c 0 >0 


Under these assumptions, H is self-adjoint in L 2 (lR n ) = Tt and its spectrum is 
pure point, and contained in ]0, oo). 

Assume that the system of non-interacting fermions is infinite and in contact 
with a reservoir at temperature T. The study of thermodynamical properties of 
this system is performed within the framework of statistical mechanics which is well 
known, and that we recall here for completeness (see [[6]]). We introduce the so-called 
Fock space: 


= ©“ „ (®:«) ( 2 . 2 ) 
where (g)™7d is the antisymmetric tensor product of Tt. which physically represents 
the space of n-fermions states. The Hamiltonian of the infinite system is governed 
by the second quantization of H : 


dT(H) = H (2.3) 

acting in T a . Similarly the number N of particles is a second quantized operator in 

T a . 


H = dT(l„) 


(2.4) 


Note that if {'0j} />o denotes an orthonormal basis of eigenfunctions of H, with 
eigenvalue Ef 
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H^j = Ej^pj (2.5) 

then: 

{Vhi A A - A ^jn} jl<j2 <...< jn (2-6) 

is an orthonormal basis of ®™7-f consisting of eigenvectors of H with eigenvalue 

Ej i + ••• + E jn- 

According to the Pauli principle, the occupation number rij of any state ifjj in T a 
equals 0 or 1. Thus the spectrum of H can be rewritten as: 

Y,j n i E j 

where rij is the eigenvalue of Nj iU : 

N jMi A ^2 A ••• A = rij (Vhi A ij) j2 A ... A i/r jn ) (2.7) 

(Nj tn “tells” whether or not the state 'ijjj is occupied in a given state of <8>"7f). We 
define: 


Nj — ® n >oEj t n 

Obviously we have: 

" = E 3 >,A- (2-8) 

H = XjEj 

Note that: [//, N] — 0 

In the grand-canonical formalism (see jSJ), the Gibbs partition function is: 

Z G = Tr (e^ S+K ^ (2.9) 

where n and f3 are Lagrange multipliers: 


f3 = 1/kT 
hi = f3(l 


( 2 . 10 ) 
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H being the chemical potential, and the Trace (which we denote with capital T) 
being taken in T a . Then it can easily be shown (see 0) that Z G factorizes as: 

= Ft , (l + (2.11) 

The mean value Fj of the occupation number of is then: 


Fj = Z G ~ l Tv (2.12) 

= (l + e^-M ))- 1 

Denoting by / the Fermi-Dirac function: 

f(x) = ( 1 + e")- 1 (2.13) 

the mean value of the number of particles in the grand-canonical ensemble is then: 

< iV >=E J > 1 tf=tr{/(/3(S- /;i ))} (2.14) 

where now the trace (which we here denote with small t) is taken in Ft. The operator 

in Ft'. 

p eq :=f((3(H-pj) (2.15) 

is called the Fermi-Dirac equilibrium one-body operator. 

We now assume that the one-body Hamiltonian is slightly perturbed in a time- 
dependent way: 


H\(t) — H + XAF(t) 

where A is a small real parameter, and F of the form, with a > 0 : 


m = { 


at 


t < o 
t > 0 . 


(2.16) 


(2.17) 


Starting at time t = —oo from the Fermi-Dirac one-body equilibrium state p eq , 
and switching on the pertubation , we get a time-dependent “density matrix” p\{t) 
(namely a trace one operator) obeying: 
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ih 


dpx 

dt 




lim p x (t) =pPg 

t^—oo 


(2.18) 


The PROBLEM is the following: to which extend does p\(t ) wander from the equi¬ 
librium state p eq as the perturbation is switched on? 


3 The linear response theory 

Physically, we aim to answer the above PROBLEM “to the first order in A”, whence 
the name “linear response theory”. In this section we give a rigorous framework to 
this program. Thus our first step is to set a convenient set of assumptions on the 
Hamiltonians under which mathematical results can be obtained. 

Assumption 2 

A(q ) is a multiplicative function dominated by 0(1 + q 2 ) in absolute value. 


Under this assumption we know (||T9|) that the unitary evolution operator gen¬ 
erated by H\(t ), namely solving: 


dVX Q t ,t,] = H x(t) Ux(m') 

v x (t 0 , to) = 1 

exists. Moreover it obeys the Duhamel’s formula: 

V\(t, t 0 ) = U(t - t 0 ) + 4 tv V x (t,lf) F(t') AU(t' - t 0 ) 

ih Jt 0 

where we have denoted: 


(3.1) 


(3.2) 


U(t) := e-" S/s . (3.3) 

We define: 

P\ (t; to) = V (t, t 0 )p eq V (to, t) (3.4) 

It is clearly a solution of (|2.18|) with p x (t 0 ) = p eq . 

We shall now justify the “linear response theory” in this context. 
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Proposition 3.1 The mapping A t—> p\(t,to ) gwen by is differentiable near 

A = 0 in the trace-class operator norm sense and we have : 

= 1 f dt' F(t') U(t-t')[p eq ,A\U(t' -t) . (3.5) 

in Jt 0 

Moreover p\(t,t 0 ) has a limit as t 0 — > — oo, in the trace-class operator norm sense , 
called p eq (t , A) and which is also differentiable in A. Moreover we have: 

= 4/ ds F (s)[p eq ,A t - s \, (3.6) 

A t being, by definition, the Heisenberg observable at time t (for the quantum evolu¬ 
tion governed by H): 

A t = U(t) A U(t)* . (3.7) 

Proof: Inserting 1 = U(t 0 — t)U(t — t 0 ) and commuting with p \/ 2 we obtain : 


d_ 

dX 


P\(t) 


d_ 

dX 


Pa (4 to 


Pa (Mo) = V x (t,t 0 )U(t 0 - t)p\ 12 .p\ ,2 U{t -t 0 )Vx{t Q ,t) (3.8) 

Each of these two factors admits a limit as to - 1 ► — oo in the norm trace sense. 
Namely using Duhamel’s formula, we have : 


(14(t, t 0 )U{t 0 - t) - 1) p\/ q 2 = A £ F (s)Vx(t, s)Ap 1 J q 2 U(s - t)ds (3.9) 

(and similarly for the adjoint) so the result follows since |P(s)|ds exists for any 
hnite t. 

Letting t 0 tend to — oo, we then get: 


^ _ A p ^ ^ 

Peq(t , A) . p eq / C?sP(s)I4(t, s) [ A , p e y]I4(s, t) 

III J — OO 


(3.10) 


A Taylor expansion near A = 0 of p eq (t, A) can be obtained by plugging in the 
Duhamel’s formula in ( [3.10D : 


V\(t, s) = U(t — s) + — I F(<r)V\(t, a)AU(a — t)da (3-11) 

ih Js 

This gives in the trace norm sense: 

Peq(t, A) - p eq = £ j dt'F(f) [p eq , A t _ t f) + o(A). (3.12) 





By this method the second term in A (“quadratic response”) could also be explicitly 
written. ■ 


Equation ( |3.10| ) is the linear response formula in this framework. It implies that if 
B is some self-adjoint operator that we want to measure in the “almost stationary” 


state p eq (t), the coefficient of the Erst order contribution in A, as A - 

0 to the 

result : 

J x (t) = tr {B ( p eq {t , A) - p e g ) 


(3.13) 

is of the form : 





J L (t ) = A f dt' F(t') $(t - t') 

J — OO 

(3.14) 

where 

= ^- tr (p [peq, A t ]) = ^tr (p eq 


(3.15) 


using the cyclicity of the trace. 


We now take the Fourier transform, in the distributional sense of <h(f), called 
the “generalized susceptiblity” : 

/ +oo 

$(f) e lwt dt (3.16) 

-OO 

which is the quantity that we shall study now. Given any function g whose Fourier 
Transform g is assumed to belong to C“(1R) : 

/ ] r+oo _ ^ ^ 

XA,B(y>)g{u)duj — — j ^ trf(f3(H - n))[A,B_ t \g(-t)dt := I(n) (3.17) 

Our aim is to obtain a semiclassical expansion of xa,b{ u) as h —> 0, f3 —> oo, namely 
a semiclassical expansion at low temperature. 

It is useful to introduce the following parameter: 

a — fdti (3.18) 

which has the dimension of time. We also define the function f a as follows: 

AM = (1 + e”)- 1 


(3.19) 
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so that I{n ) can be rewritten formally as: 

m = W, / tr {~h^) lA ’ St] ) mit (3 ' 20) 

However this expression suffers from the singularity in 0 of f a as a —> 0. In order 
to avoid this, we “regularize” it by using instead of f a \ 


fa, V = fa*V (3-21) 

where r] is a function in 5(R) such that its Fourier Transform rj G Cq°(R). This 
amounts to study Xa,b as a distribution on R 2 in the variables s and u in the 
following way: 


1 

2i7rh 


j j XA,B{s,uj)r](s)g(uj)dsduj = 


tr f e is{H ~^ /n 


A,B t ) U(s)r](s)g(t)dsdt 


(3.22) 


Let us intoduce the following test space functions on R 2 = R s x R w : 

Definition 3.2 We say that p G K a , a > 0, if p is C°° on R 2 and there exist b > 0, 
c > 0 such that p(s,u) = 0 for |s| > b, uj G R, and |<^ 2 )(s, t)| < ce~ a ^ for every 
(s,t) G R 2 , where <p^(s,t) denotes the Fourier transform in the second argument. 

4 The results 

In this section we first introduce the notations of the classical objects that will 
appear in the semiclassical expansions, together with the assumptions under which 
these expansions can be obtained. 

Let (ff be the classical flow induced by Hamiltonian 
surface conserved by the flow: 

^ = {(q,p)eB 2n :H(q,p) = p} ( 4 . 1 ) 

We call dEp the Liouvillc measure on E ;i , so that the correlation of classical observ¬ 
ables A and B on E^ is defined by: 

CA,B,n(t ) = I A.B t dY>n 


(2.1). Consider E M the energy 


(4.2) 
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where B t {z) = B[<p t {z)]. Moreover if 7 is any periodic orbit on , and 7 * the 
corresponding primitive orbit, with period T 7 *, we introduce the correlation function 

r'l 7 * 

Cj*[t)= A s (q,p)B s+t (q,p)ds (q,p) e 7 * (4-3) 

Jo 

c 7 * being T 7 *-periodic it admits the Fourier-series expansion: 

k=+oo 

c 7 *(t) = c^ yk e 2lktn/Tnl * (4.4) 

k=—oo 

To each 7 is associated a corresponding “linearized Poincare map” called P 7 , a 
classical action along 7 called S 7 , and a Maslov index z / 7 (see 0 ). 

Let us assume that (j) t on satisfies the so-called Gutzwiller Assumption: 
Assumption 2 The periodic orbits 7 are non-degenerate, i.e the Poincare maps do 

not have 1 as eigenvalue (which implies that they are isolated). 

Moreover we shall be able to treat B obeying: 

Assumption 3 

\^a^B(q,p)\<C a g |a| + |/3|>2 

Our result is as follows: 

Theorem 4.1 Under Assumptions 1,2,3, we have, in distributional sense in JC a 
(see definition \3.2p : 

Xa,b{s,u) = - h~ n 5 0 (s ) < 8 ) C' AtBtli (u) + ^2tv , ~ n fi j (s,u) 

j> 1 


+ E 

r-T-y^o 


7t/2) 

ha sinh (nT^/a) \ det(l — P 7 ) I 1 / 2 


^T-y (S) 


k 




6(u 



+ E^E>7(A^) 

j> 1 


+0{h aiH - £ ~ n ) 

where pj and u, jri are distributions in JC a such that Supp(pj) C {0 }xR, Supplyj^) C 
x It, and 7 h is a non negative constant depending only on H and p (not on 
a). 


Proof: 

As a distribution acting on (rj ® g)(s,u), Xa,b is given by (|3.22| ). We split the 
integral over t into two parts: \t\ < 7 nLog(l/h ) and its complement, where 7 h 
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is a constant obtained in Egorov-type estimates (see f|) and only depending on 
Hamiltonian H. 


Using the exponential decrease of g(t), it is not difficult to estimate the contribution 


of the integration domain \t\ > 7 nLogil/h ) as 0(h a " /H £ n ), for any e > 0. The 
larger is a (the exponential fall-off rate of g) the smaller is this “error term”. 


In order to estimate the contribution of the integration domain \t\ < 7 nLog{\jTi ) 


we shall use truncations in the spectral variable of Hamiltonian H in order to apply 
known results and usual methods. 

In all that follows, the integration support in t variable is supposed to be \t\ < 
'ynLog(l/h ), and we call I v (g) the resulting contribution to ( |3.22| ). Fix 5 positive 
and small enough and let us introduce a C°° partition of unity as follows : 


1 — C- + Co + c+ 


(4.5) 


where 




(4.6) 


and SuppC- C] — ex), —5/2], SuppC+ Q [5/2, +oo[. 


Inserting in (|3.22|) 


l = C-(H-g) + (o(H ~n) + C+(H - n) 


we obtain, with obvious notations : 


W = W +W + VM 


(4.7) 


Let 9 be a regular Schwartz function such that its Fourier Transform 6 be in C/°(1R,), 
and 



if \t\ < 1 
if \t\ > 2 


(4.8) 


For any positive number r, we set : 


0 T (s) — 8 {s/t) . 


(4.9) 
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and let us denote by 0 T the inverse Fourier transform of 0 T . Let To be a positive 
number, small enough in a sense to be made precise later. We shall now decompose 
7°(/i) in two parts: 


W = + 


2iTrh 


tr 


J j dsdtf a (s)T](s) (l - 0(s/ T 0 )) e" {H “ )/h [A B t ] Co (H - ft)g(t) 


(4.10) 


Thus equ. (|4.7|) now becomes: 


40) = CO) + CM + 4 + M + 4" (a) (4.ii) 

where each term can be estimated separately. 

Estimate of : 

Denote by : 

0/3® :=/(/?(£-/r))C_(£-h) 

We remark that 0oo = C~(E — //), and 

0/3® = (4.12) 

for some x G C^°(1R) because the spectrum of H is bounded from below. 

But E i—>• (j)p(E)x(E ) is a bounded family of functions in Qj°(ZR) for f3 in ]0, +oo]. 
Therefore the U -semiclassical functional calculus can be applied, yielding an asymp¬ 
totic expansion of the following form : 


4“M~R‘”Ecj^ 

j>o 


uniformly in a e]0, +oo], where : 


(4.13) 


c 0 = ( 2 tt) ” j f ifi{H(q,p) - n))(-(H(q,p) - n) {A,B t } (q,p)dqdp (4.14) 

and analogous formulae for j > 1 , where we have used the known result that the 
principal symbol of: 


Ct 


i 

h 


[A St] 


(4.15) 
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is {A,B t }, using Egorov’s theorem in the form given in @4[] for \t\ < ^nLog^l/h), 
and semiclassical calculus. 

Estimate of /+ (/i) : 


IC»I < CY. I fWE, - »))U E, - n)\ (4.16) 

... 

where C depends on g and A, B , and where {Ej} . >± is the increasing sequence of 
the eigenvalues of H 


<C E l/m-K))l (4.17) 

Ej>5/2 

We introduce the counting function : 

N(E ): 

Using a “Lieb-Thirring-like” estimate, we get : 

N(E) < r yh~ n (l + E) m 

We therefore deduce the existence of a positive constant c (depending on 5) such 
that : 


ICJ < Ch -"e- c|3 (4,18) 

for any h e]0, 1] and any positive a. 

Estimate of /“ TQ (/i) 

Recall that : 




Vo(^) = 1 / ih i /_ oo tr ^,*7, 

where we have defined: 


TO 


A, B, 


g(t)dt (4.19) 


fcr,ri,To ft t * T) * 9 


ro 


(4.20) 


We want to estimate 


L,y-= ftr ( 0 (H - ft) 


A.B, 


(4.21) 








using equ. (|4.11|) we can write : 



Co(H - ii)C t d\ 
(—^—) Co(H — n)(3 t d\ 


(4.23) 


We have thus: L Mjt = L 1 1 + L 2 ^ t , with 


\Llt\<Ch~ x 


rii—25 


E 

3> 1 


(K*v*Qt 0 )( Ej A 


h 


)Co (Ej - fj) 


dX 


(4.24) 


where C is uniform with respect to t G Suppg and to h. 

By playing with localization and decay properties, one easily obtains that L 2 t = 
0(h°° ) uniformly with respect to t G Suppg, and with respect to a e]0, +cxd[. 

The term L^ t can be dealt with as in ||, using the fact that for 6 small enough, A 
is non critical for H for every A G [/x — 26, /i]. Thus L l t can be rewritten as: 


L X = -f JW iX 

J \l— 2d 

where : 


(4.25) 


I(\)=ih 1 / is TA 9{s/T 0 )ri(s) tr(e (4.26) 

J-oo smh S7T/(T l J 

Then using a coherent states decomposition of the trace as in [|7|, we see, using 
the support property of 9 ro that the dominant contribution in the stationary phase 
theorem comes from s = 0. Thus, provided that 2tq is smaller than the smallest 
period of closed orbits on equ. ( |4.26| ) provides an asymptotic expansion in h, of 
the form: 
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/(A) = h~ n (C 0 ( A) + hCi (A) + ....) mod 0(/i°°) (4.27) 

which is uniform in cr e]0,+oo[, and which can be further integrated with respect 
to A on the interval \p, p + 25], yielding the result. 


We shall now give the explicit form of the dominant 0{Ti ~ n ) contribution to /^(/i) 
(equ. ( |4.6|) ) which comes from the sum of the contributions of J~(p) and /” To (/z); 
we obtain: 


= h n r](0) [ g(t) [ {A, B t } (q,p)dpdqdt + O^h 1 n ) (4.28) 

We have introduced the correlation in time of A and B on the energy surface 
= [H(q,p) = n\ (see flU)) 

-■=L, AB ‘m (4 - 29) 

Let be a C°° function with compact support contained in ] — oo, p + 5]. We have : 


J{A,B t }p(H(q,p))dpdq = J {Ap(H), B t }dqdp - J A{p(H), B t }dqdp (4.30) 

where the integration is over the full phase space !R 2n . Bp(H) being a Q(°(R 2ti ) 
function of ( q,p ), we get by integration by part that : 

J{A<p(H), B t }dqdp = 0 

Moreover 


J{Ap(H), B t }dqdp = J {H, B t }<p'(H)Adqdp = J AB t tp\H)dqdp 
We now let <p tend to and find : 

[ {A, B t }dqdp = ( f A.B t -^j—\ 

Therefore the dominant term of J r; (/z) is given by : 

IM = h~ n m J C A ^(t)g \t)dt + 0(h l ~ n ) 
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This completes the proof for the first term of the asymptotic expansion in The¬ 
orem (4.1). 

Estimate of I° s f TQ (/i) 

We have : 


/“» = h- 1 j digit) r{Co(if-M)e-" (S -" )/, C,} (4,31) 
where we have used the following notation : 

Ve,T 0 (s) ■=rj(s ) (l - Ms)) 

Again we proceed as in [7] by a “Gutzwiller type” estimate for the integral over s 
since the support of Tj gT0 (s) doesn’t contain s = 0 , but will only contribute by a 
finite number of closed classical orbits which we denote by 7 . Furthermore due to 
the support properties of 6, it is clear that r]e tT0 (T^) = 7 ](T 1 ) 

Using the Gutzwiller assumption, and the compact integration support in 
variable t, we obtain the following asymptotic expansion of /"^(/i), which is uniform 
in the parameter a e] 0 , + 00 [: 


rose 

1 rj 


(to (aO = h 1 1 dtg(t) J 2 


7 r 


giS~f /h+ii/'y 7 t /2 


0 sinh 7r Tj /<7 | det(l — P^l 1 / 2 


x 


(4.32) 


x ( 'Ev(T 1 )c r , k e 2i * kt/T l + 5 >i/ jVy (t) 

; j>i 


This yields the following contribution of oscillating terms to xa,b(s,u>), in the dis¬ 
tribution sense, uniformly in a e]0, +00 [ : 


ire 


i(S-y/h+u^ir/2) 


Xa,b,os C ftcrsinh (7iT^/a)\det(l — P^l 1 / 2 

/ _2/C7T _> 

x ( £t 7 (s) ® ^c 7 *, fc h( U- — ) + J2 hJjy j,'y(s,uj) 


' y* 


j>i 


(4.33) 


where Uj tJ (s,u) are distributions supported in {T 7 } x 1R 
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5 Concluding remarks 


Theorem 4.1 is an extension of the well known Gutzwiller trace formulae for the 
spectral density of energy levels. The main difference is that here there are two real 
variables instead of one because in the “dynamical susceptibility” time and energy 
variables are mixed up in an intricated way. So we can put in a mathematical 


rigorous shape the main result of the paper J13 


As in the Gutzwiller trace formulae our formulae in Theorem 4.1 gives a semiclassical 
expansion with three different terms: the first line gives a regular expansion in h, 
which is the contribution of the period 0 of the classical flow; the second line is an 
oscillating part coming from the contributions of the non zero periods of the classical 
flow; the third line is the error term depending on the test functions considered. 

So far we have shown that a semiclassical expansion, in the linear response the¬ 
ory, can be obtained for a regularized version of the “dynamical susceptibility”, i.e 
in a suitable distributional sense. The same is obviously true for the linear response 
function Jhif) defined by ( [3.15|) , as we shall establish now. 


Formally, if tp is a C“(R) function, we have, in distributional sense: 

< J L ,tp >= f h(u })xA,B( U ) duJ 

where ki(u) := Q(u)k(u), 

(0 being the Heavyside function) 
and k(u) := f p(s + u)F(s)ds 

However Xa,b is only well defined mathematically as a semiclassical expansion in a 
“regularized” form: 


Xa,b, v ■= f r](s)xA,B(s,uj)ds 

where i) is in C^°(R). 

Similarly, a “regularized form” of Jl can be defined as: 

JL, v (t) := J fj(s)J L (t,s)ds 

in the following sense: 




ki(u)xA,B,r 1 (s,u)du 
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So we have 

< J L ^ v (f >= J fj(s)k 1 (u)xA,B( s ^ u )dsdw 

It is not hard to see, using the definition of k\ that rj(s)ki(u) G /C Q _ e for any £ > 0, 

so that our theorem applies. For example we can compute the leading term : 

< Jl n-, T >= — 2nh~ n fi(0) I duC' AB (u ) [ dsip(s + u)F(s) + 0{h l ~ n ) 
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